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Abstract 



Let a "complex probability" be a normalizable complex distribution P(x) 
defined on K D . A real and positive probability distribution p(z), denned 
on the complex plane C D , is said to be a positive representation of P(x) if 
(Q(x))p = (Q(z)) p , where Q(x) is any polynomial in R D and Q(z) its analyti- 
cal extension on C D . In this paper it is shown that every complex probability 
admits a real representation and a constructive method is given. Among 
other results, explicit positive representations, in any number of dimensions, 
are given for any complex distribution of the form Gaussian times polynomial, 
for any complex distributions with support at one point and for any periodic 
Gaussian times polynomial. 
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I. INTRODUCTION 



In quantum physics there are instances of averages where the role of probability distri- 
bution is played by a distribution taking complex values. Consider the functional integral 
formulation of field theory There, the time ordered expectation value of observables 
takes de form (TO[(j)]) = N J T)(p(x) e lS ^0[(p], where S[4>] is the action functional and iV 
a normalization constant. This is a first instance of a "complex probability distribution", 
namely, the Boltzmann weight P[(f>] = Ne lS ^. In the continuum, such functional integral is 
not sufficiently well-behaved and only its Euclidean version can be given a rigorous meaning 
0. Within a lattice regularization, the Minkowski version is mathematically well-defined, 
nevertheless the Wick rotation is performed in this case too. This is because, in most cases, 
in the Euclidean theory the Boltzmann weight becomes a real and positive probability dis- 
tribution. This is important in practice since straightforward Monte Carlo is only defined for 
positive probabilities. There are cases, however, when even Euclidean field theory presents 
complex actions. Indeed, the statistical interpretation of the quantum theory requires the 
Boltzmann weight to be reflection positive, but not directly positive ||. Instances of com- 
plex Euclidean actions occur after integration of fermions, since the fermionic determinant 
is not positive definite; if there are non vanishing chemical potentials; in gauge theories in 
the presence of Wilson loops or topological ^-terms or in general after inserting projection 
operators in the path integral to select particular sectors of the theory |^,[5,[lT,7|. Also, two 
dimensional fermions can be brought to a bosonic complex action form ||. 

As we have said, the computation of averages in the presence of a complex probability 
distribution poses a practical problem, namely, the Monte Carlo method cannot be used 
directly to sample the probability since this method only makes sense for true, i.e. real and 
positive, probabilities. The standard approach to complex probabilities in numerical simu- 
lations [|],|5| is to factorize a real and positive part to be used as input for some Monte Carlo 
method and include the remainder in the observable. That is, if the complex probability is 
P(x) = Pq(x)F(x) with Po(x) positive, the expectation values can be obtained as 

Of course, the same formula can be used when P(x) itself is positive. The problem with 
this approach is that it violates the importance sample principle, since we are not sampling 
the true probability and that increases the dispersion of Monte Carlo data. For instance, 
(F(x))p Q may be small, thereby introducing large error bars. 

An alternative approach is to look for a positive probability p(z) in the complex con- 
figuration space which gives the same expectation values as P(x), i.e., {0(x))p = (0(z)) p , 
where O(z) is the analytical extension of 0(x). The usual way of constructing such a 
probability is by means of the complex Langevin algorithm P,10|. In this approach the 



configuration is updated through a standard Langevin algorithm with the complex action. 
Since the drift term is complex, the complex extension of the configuration space is sampled 
as well. Whenever the random walk possesses an equilibrium configuration, it is sampling 
the complex configuration space with a real and positive probability distribution p(z). We 
have then traded a complex probability P{x) on R D by a positive probability p(z) on C D . 
If p(z) happens to be equivalent to P(x) in the sense of expectation values, we have suc- 
ceeded in sampling the complex probability. Successful implementations of the algorithm 
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have been obtained in some practical cases, such as two dimensional compact QED with 
static charges ||. In general, however, the complex Langevin algorithm poses two problems. 
First, it not always converges to an equilibrium distribution. Second and more subtle, for 
some actions it seems to converge to an equilibrium distribution which is not equivalent to 
the original complex probability [|TT|-|T3|j, (see however [Q). Such phenomenon has been 
found in practically relevant cases such as QCD with a Wilson loop [|TT| , |l^ .|T5|] . 

In the present paper we consider the problem of constructing a positive representation 
directly, independently of the Langevin algorithm. Several properties of representations 
of complex probabilities on R D by probabilities on C D are noted. A constructive method 
is given to obtain real (although not necessarily positive) representations of very general 
complex probabilities. Positive representations are explicitly constructed for some probabil- 
ities which are beyond the present applicability of the complex Langevin algorithm. These 
include Gaussian times polynomial, distributions with support at one point, and periodic 
Gaussian times polynomial. In all cases, such representations are not unique. 

These results are of great interest from the point of view of applications. This is not 
because the constructions found here are of direct usefulness to carry out numerical cal- 
culations; there are far more natural ways to compute expectations values with complex 
Gaussian times polynomial distributions. The interest lies in the following. The negative 
results found up to now with the complex Langevin algorithm in some systems would make 
one to have reasonable doubts of whether a positive representation exists at all for those 
systems. Moreover, the momenta of any positive probability on C D are bounded to satisfy 
some inequalities among them. It might happen that those bounds were incompatible with 
the momenta of the given complex probability on R D in some cases. At present, the nec- 
essary and sufficient conditions for a positive representation to exist are not known. The 
results of this paper suggest, however, that such representation exists quite generally since 
the set of Gaussian times polynomial is dense in L 2 (R D ). Our results tend to support the 
idea that there is no obstruction of principle for positive representations to exits. This is 
the main insight of this work. 



II. REPRESENTATION OF COMPLEX PROBABILITIES 

The complex probabilities P(x) to be considered here will be tempered distributions 
on K D of a restricted class, namely, those which are the inverse Fourier transform of an 
ordinary function P(k) (locally integrable and at most of polynomial growth at infinity), 
with P(k) non vanishing at the origin and analytical at that point. These conditions allow 
for a natural definition of J x^ ■ ■ ■ Xi n P(x)d D x through the Taylor expansion of P(k) at 
k = 0. In particular / P(x)d D x will be non vanishing. The expectation value associated to 
P(x) is defined for any polynomial Q(x) as 



jQ(x)P(x)d D x 
JP(x)d D x 



(Q{*))p = J T Z ,0X ■ ( 2 ) 



Likewise, we can consider complex probabilities on C D as the class of distributions defined 
above on R 2D . For any such distribution, p(z), the expectation value takes the form 
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p 



fq(z)p(z)d z 
Jp(z)d 2D z 



(3) 



where Zj = Xj+iyj, d z = d xd y and q(z) is an arbitrary polynomial of z and its complex 
conjugate z*. 

By definition, p(z) is a representation of P(x) if (Q(x))p = (Q(z)) p , where Q(x) is 
any polynomial on K D and its analytical extension on C D . Equivalently, one can 

demand (x^ ■ ■ ■ x in )p = {z ix ■ ■ ■ z in ) p for any set of indices, where i r = 1, . . . ,D and n = 
0, 1, 2, . . .. Two complex probabilities on C D will be called equivalent if they have the same 
expectation values on every analytical polynomial. In general, two equivalent probabilities 
will not coincide on expectation values of non analytical polynomials (z^ ■ ■ ■ z in z* x ■ ■ ■ z* m ). 
A representation will be called real if p(z) is real, positive if p(z) is non negative and unitary 
if J p(z)d 2D z = J P(x)d D x. Our goal is then to find positive representations of complex 
probabilities. 

We will proceed by noting different ways to obtain new representations from known ones. 
A first obvious way is by means of complex affine transformations. Let A be a non singular 
complex D x D matrix, and a £ C D , and assume that Pq(z) is an analytical function 
in a region including R D and AR D + a such that Po(x) and P(x) = det(A)P (Ax + a) 
are both complex probabilities. Then if po(z) is a unitary representation of Po(x) so is 
p(z) = | det(A)\ 2 p (Az + a) of P(x): for any polynomial Q(x) 



Furthermore, p(z) is positive if po(z) is positive. Another construction follows from linear 
combination. If Pi{z) are unitary representations of P%{x\ so is p{z) = Y^i=ihPi(z) of 
P(x) = Yh=i biPi(%)- Again, if Pi(z) are positive and bi non negative, p(z) is positive too. 

Let us define the partial derivatives dk and &l on a function on C D as (d / 'dx^Tid / 'dyu) /2, 
respectively and let 4>{z) be in the class of distributions on C D defined above but dropping 
the restriction / (f)(z)d 2D z ^ 0. Then if p(z) is a probability, p(z)+dl<j)(z) is also a probability 
and in fact (unitarily) equivalent to p(z), 



where Q(z) is any analytical polynomial. That is, dl<p(z) would represent the zero distri- 
bution on K D . Such distributions will be called null distributions. They will prove useful 
in what follows to obtain positive representations from real ones, namely, by adding null 
distribution of the form J2k=i dk&k^kiz) ■> f° r suitably chosen real <j>k{z). Note that 4<9fcC^ is 
just a Laplacian. 

Similarly, by proceeding as in eq. @, it follows that if p(z) represents P(x), the following 
relations hold 



det(A)\ 2 J Q(z)p (Az + a)d 2D z = J Q{A- 1 {z - a))p (z)d 2D z 

= ( Q(A~ 1 (x - a))P (x)d D x = det(A) [ Q(x)P (Ax + a)d D x 



(4) 



J Q(z)d* k( l ) (z)d 2D z = I dl{Q{z)^{z))d 2D z = 0, 



(5) 




(6) 
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where Q(x) and R(x) are arbitrary polynomials. That is, dk on C D represents dk on R D and 
multiplication by an analytical polynomial R(z) represents multiplication by R(x). 

Another interesting construction is related to convolutions. The convolution exist for 
any two complex probabilities since it can be defined through the product of their Fourier 
transforms which are regular distributions. If pi(z) and p 2 (z) are unitary representations of 
P\{x) and P2{x) respectively, their convolution pi *p 2 is a unitary representation of Pi * P 2 . 
Indeed, p\ <8> P2 is a unitary representation of Pi <S> Pi and 

( z i\ ' ' ' z i„)pi*p2 = (( Z iy Z ii ) ' ' ' ( Z i n Z i n ))pi®P2 

= + • • • ( X S + 4! } ))Pl®P 2 = fail ■ ■ ■ Xin)Pl*P 2 ■ (7) 

Furthermore, if p±(z) and p 2 (z) are positive, pi *p 2 is positive too. In particular, this allows 
for obtaining equivalent representations of known ones: if p(z) is a unitary representation of 
P(x) and C(z) is a unitary representation of 5(x), the .D-dimensional Dirac delta function, 
p*C will be unitarily equivalent to p(z), since P * 5 = P. Any probability C(z) normalized 
to one defines a unitary representation of 5(x) if it is invariant under global phase rotations, 
i.e., C(e ltp z) = C(z) for any ip e R. In this case 

]2D 



h ---z in C(z)d 2U z = 5 n , , (8) 

since the angular average of Zi x • • • Zi n vanishes for n > 0. In fact this construction can be 
regarded as adding a Laplacian, namely, p*C — p, as it is easily seen after Fourier transform. 
This procedure can be used to obtain positive representations from real ones. On the other 
hand, it shows that if a complex probability admits a unitary positive representation it is 
not unique. 

A unitary representation can always be obtained for any P(x) by taking p(z) = P(x)5(y). 
If P(x) is positive so will be p(z). This can be generalized as follows. Let Po(x) be positive 
and P(x) = Pq(x — it), t e R D (i.e., a complex translation under the conditions considered 
above for affine transformations). Then p(z) = Po(x)5(y — t) is a unitary positive represen- 
tation of P(x). If we allow P to depend on t, taking linear combinations we obtain that 
p(z) = p(x, y) is a unitary representation of 

P(x) = J p(x - iy,y)d D y . (9) 



This relation has been noted before in the literature [T^|T^|, considered as a projection from 
probabilities on C D to probabilities on R D . Note, however, that when this relation can be 
applied it gives just one of the P(x) represented by p(z). In fact, since the momenta of 
P{x) are the Taylor expansion coefficients of its Fourier transform, there are many complex 
probabilities characterized by the same momenta. As we have seen, under this projection, 
the operation d* is mapped to zero. Similarly, <9j is mapped to d/dxi, and multiplication by 
an analytical polynomial Q(z) is mapped to multiplication by Q(x). 

As an immediate application of eq. (|^), we find that for rriy real, symmetric and positive 
definite, the probability 

P(x) = /(x)exp(-^m ij x i x i ) (10) 
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is represented by 



p(z) = det(m)f(my) exp^m^z*) , (11) 



where / is the Fourier transform of / (the repeated index convention will be used in what 
follows). For example, for D — 1, and T positive, P(x) = cos(a;) exp(— x 2 /2T) is represented 
by the positive probability p(z) = exp(—x 2 /2T)(5(y — T) + 5(y + T)). Since in this example 
P(x) is real but not positive definite, this is an instance where a complex Langevin simulation 



would fail |T^ , p!5|JT3| , yet there is a positive representation. 

Next, let us show that every complex probability on K D admits a real representation. 
Let P(x) be a complex probability normalized to one and P(k) its Fourier transform 

p(k) = J e ikx P{x)d D x (12) 

where kx = kiXi. By definition we have 

°° i n 

?( k ) = Yl —\ k n ■■■K (x h •■■x in ) P (13) 



n=0 U - 



in a neighborhood of k = since P{k) is analytic at the origin. Also, 

(x h ---x in ) P = (-i) n d h ---d in P(k)\ k=0 . (14) 
For a probability p(z) on C D , the Fourier transform is defined similarly, 

p(a) = J e ikx+iry p{z)d 2D z = J e i{a * z+az * )/2 p{z)d 2D z, (15) 

where <ii = ki + iri. Assuming that p(z) is normalized to one, its momenta are obtained 
through 

(*i ■ • • ^4 ■■■ z Vp = i-^T^dl ■ ■ ■ d; n d n ■ ■ ■ d lm p(a) \ a=0 (16) 
where <9j refers to ex, and d* to a*. Consider the following probability, 

p{a) = C{a)pC-){P{- G -)r. (17) 

Here C(z) is one of the real unitary representations of S(x) above mentioned. Thus C(cr) 
is analytical at the origin as a function of ki and rj and is invariant under global phase 
rotations of a. P(<j) stands for the analytical extension of P(k) in a neighborhood of the 
origin. Beyond the analyticity circle (if it is finite) we can choose C(a) equal to zero so 
that p(a) exists. By construction, p(a) is unity at the origin and analytical there. Also it 
is locally integrable and, with a suitable choice of C{<j), grows at most polynomically at 
infinity, therefore it defines a probability p{z) on C D . Furthermore, p(z) is real since C(z) 
is real and (p(cr))* = p(—a). It remains to show that it is a representation of P(x), 



• z in ) P = (-2i)X " " " d;j(a) 1 =q = (-2i)X " " " ^P(y) 



<T=0 



<-i)"d h ---d in P(k) k=o =(x il ---x in ) P , (18) 
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where it has been used that d* ■ ■ ■ d* n C(a)\ a=Q vanishes for n > 0. That is, we have given 
a constructive method, eq. (|17D , to obtain a real representation of any complex probability 
within the class of complex probabilities considered. 
As an illustration, consider D = 1 and 

P(x) = 5(x) + a5'(x) , a = clr + iaj £ C . (19) 

In this case -P(o~) = 1 — iaa is a polynomial, thus it is entire and well-behaved at infinity 
and we can take C{a) = 1, i.e., C(z) = S(x)S(y). With this choice 

1 i 

p(a) = 1 - -|a| 2 |a| 2 - -(aV + acx*) (20) 

and 

p{z) = 5(x)5(y) + a R 5'(x)5(y) + aj5(x)5'(y) + ]\a\\5" {x)5{y) + 5(x)5"(y)) . (21) 

One can easily check that this is a real distribution which represents P(x), however it is 
not positive. We can find a positive representation by first applying a convolution (i.e., a 
better choice of C(z)) and then adding a suitable Laplacian. Furthermore, it can be done 
for an arbitrary distribution of support at zero in any number of dimensions. Rather than 
showing this in detail here, it will be obtained as a byproduct in the next section. There we 
will obtain positive representations of Gaussian functions times polynomials. 

By formally undoing the Fourier transform of p(a) in eq. fllTD , the following explicit form 
of p(z) is obtained 

p(z) = J C (x - ^±El, y - ^-^)P( Xl )P*( X2 )d D x 1 d D x 2 , (22) 

where C (zi, z 2 ) is the analytical extension of C (x,y) = C(x + iy), with x and y real. In 
order for this formula to make sense, we should require Cq(zi, z 2 ) to be entire on C 2D and 
further the integrand should be sufficiently convergent so as to define a probability on C D . 
Such probability is real by construction, since C(z) is real, however it will not be positive 
in general even if C{z) is positive since such property is lost after analytical extension. 
The interest of this relation, as compared, for instance with that in eq. (Q), is that it is 
constructive. 

An example of application of this formula is provided by 



N 



P(x) = £ ai 5{x - x«), C(z) = exp(-^f ) , (23) 



i=l 

which gives 

p(z) = a ia * exp f -- f \x J +[y ) | | . (24) 

Another application is when P(x) is a finite linear combination of Gaussian distributions 
centered anywhere in the complex plane and with arbitrary complex widths, provided we 
choose r > \Ti\, i — 1, . . . , N . 
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III. POSITIVE REPRESENTATIONS OF GAUSSIAN DISTRIBUTIONS 



A Gaussian complex probability takes the general form 

G(x) = Ncexpi-^rriijXiXj - 6;X;) , 
N G = (27r)- D / 2 (det(m)) 1 / 2 exp(-^(m- 1 )^6 J 6,) (25) 

where symmetric complex matrix with positive definite real part to ensure normal- 

izability. As a consequence is non singular and can be written as A^A^. This allows 
to set rriij = 5ij and hi = by means of a complex affine transformation. That is, we will 
consider only 

G(x) = (2n)- D / 2 exp(-^x i x l ) (26) 

and the general case can be obtained a posteriori as G(Ax + A^ 1 b). A positive representation 
of G(x) is simply G(x)5(y). A more general representation g(z) is obtained by convolution 
with C(z) = (27n])~ D exp(— ^ZiZ*), where 77 is positive. This gives 

1 1 
9\ z ) = ^exp(— — — — XiXi - —ViVi) 
2(r] + 1) 2r] 

1 1 

= N g exp(-7Zi< + -jZiZi + -jz*z*) , (27) 

where the normalization constant is N g = (271-^/77(77 + 1)^ and we have introduced the 
positive numbers 

1 277 +1 

7 = , 7 = — • (28) 

' 477(7/ + 1) ' 477(7/ + 1) ( J 

The same representation is obtained by following the method of eq. (^). The value of the 
parameter rj, or equivalently 7, will be fixed below. 

The set of probabilities to be considered is P{x) = Q(x)G(x), where Q(x) is a complex 
polynomial of degree N. P(x) can always be written as 

N l 

P(x) = 52-a ll ^ n a il ...a in G(x), (29) 

n=0 n - 

where a^...^ is completely symmetric and the zeroth order coefficient a must not vanish (in 
fact, is unity if P(x) is normalized). A real representation of P(x) is given by 

N -y N 1 

Po{z) = (\a Q \ 2 + a*Y, -ja h ... tn d h ■ ■ ■ 9^ + a ^ —a* ... in <9* . . . d*Jg(z) , (30) 

n=l n - n=l n - 

since the terms with d* do not contribute and d/dz is mapped to d/dx under projection. 
It is convenient to introduce the polynomials 
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Qh...i n (z) =g(z) 1 d h ---d in g{z). (31) 

They can be computed recursively by means of the formula 

Q (z) = 1, Q h ...i n ( z ) = (d in + Wi n )Qi 1 ...i n - 1 (z) , (32) 

where we have introduced the variable 

Ui = - "yz* . (33) 

The functions Qi 1 ...i n (z) are polynomials of degree n in u^, with coefficients depending only 
on 7. With this notation, po(z) can be rewritten as 

/ N 1 N 1 \ 

Po(z) = \a \ 2 + a* ]T — ] a h ...i n Qi 1 ...i n ( z ) + a o E Z^L.^L.^) ^ ■ ( 34 ) 

V 71=1 71=1 / 

In order to obtain a positive representation, p(z) can be further cast in the form 

p (z) = ( I « 1 2 + E —Pn\Qh...i n (z) + (3n la 0a* h ...iJ 2 



n=l 
N 



-E(^l^l...^W| 2 +/3- 1 |«0| 2 |«n| 2 ))^W, (35) 



n=l n '- 



where the indices i\ . . . i n are summed over. f3i, . . . , Pn are arbitrary positive numbers which 
value is to be specified below and we have defined the quantity | (Xyi I clS 



,2 1 



Am. ..>,,< ...,„• (36) 



We will assume that \a n \ is non vanishing, since the vanishing case is trivial. In Appendix 
A it is shown that 

Uz) = (^Qn..,MQl..,S z ) - f l K n (D^ g{z) (37) 
is a null distribution, where 

Kn{D) ~ n\{D-l)\ ■ (38) 
By removing (f) n (z) from po(z) we obtain an equivalent representation p(z), namely, 



p(z) 



JV 1 N 

E -Pn\Q ll .., n (z)+fa 1 aoat..J 2 + \a \ 2 - E^f 1 *^) + (3- l \a \ 2 \a n \ 2 ) 

.n=l U - n=l 



</(*)■ 

(39) 

To ensure positivity of p(z) we require 



TV 



Y(Pnl n K n (D) + ^-VoHanl 2 ) < |a | 2 . (40) 



71=1 
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This can be achieved by choosing the positive coefficients (3 n so as to minimize the left-hand 
side, 

Pn= ' a ° IH ■ (41) 

H n K n (D) 



In this way the inequality is satisfied for any 7 smaller than the unique positive solution of 

N 1 1 

J2^K n (D)\a n \r /2 = -\a \. (42) 

71=1 Z 

For this choice of 7, p(z) takes the simple form 

N j 

p(z) = E -Pn\Qn... in (z) + f3- l a Q a*^ i fg{z) . (43) 

71=1 U - 

To summarize, any Gaussian times polynomial complex probability, eq. (p9p, admits a pos- 
itive representation, namely, in eq. (f4"3"D, with /3 n given by eq. (0), and 7 given by 
eq. (H). 

Incidentally, let us note that from a computational point of view, it is convenient to 
minimize the width of p(z) in the complex plane (e.g., if P(x) is already positive, the 
best choice is P(x)5(y)), since this reduces the dispersion of points in the sample. In the 
family of probabilities described by the expression of p(z) in eq. fl39|), this minimization 
corresponds to our choice of f3 n in eq. (|4l|) and 7 in eq. (f42"l). In general, however, this needs 
not be best equivalent positive representation of P(x). The construction presented above 
corresponds to adding to po(z) a Laplacian of the form ■ ■ -d^d^ ■ ■ -d* n g(z) (as can be 
seen using the formulas of Appendix A). More generally, one could add terms of the form 
h 1 ..A„;j 1 ...j n di 1 ■ ■ ■ d^d^ ■ ■ ■ d* n g(z), with b self-adjoint, in order to optimize p(z), or even more 
general terms so long as they have a d* and are real. 

Let us now come back to the problem of finding positive representations of complex 
distributions with support at 0. Such distributions take the form 

N 1 

P( x ) = J2^a h ... in d h ...d l J(x). (44) 

This distribution can be considered as the zero width limit of the Gaussian times polynomial 
distribution. 

N 1 

P(x) = \imP x (x), P x (x) = Y J -^..,A,---d in {^ D G{x/\)). (45) 

Naming P(x;a) the probability in eq. (P9"D, we find 

P x (x) = \- D P(x/X; a x ) , al_ An = A~X...*„ • (46) 

Therefore, the positive representation of P(x;a), namely, p{z\ a) in eq. ([43D , provides a 
positive representation of P\(x), 
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p x {z)=\- w p{z/\ ] a"). (47) 

In order to take the limit, we should consider how the different variables scale. We already 
have the scaling law of z and of the coefficients a^...^. From eqs. (|4l|j4^) /3„ is found 



to scale as \~ 2n f3 n and 7 A as A 2 7- From eqs. (|28|), r] x is given in leading order by \~ 2 r] 
with r] = 1/(27) an d 7 A is °f order A 4 and can be neglected. Therefore, in leading order 
\~ 2D g(z/ \; 7 A ) becomes 

g\z- 7) = (27r V )- D exp(-7^<) , 77 = (48) 

27 

and is independent of A. This results is to be used in eq. (^). Finally, in leading order, 
Qn...i n {z/\;i X ) becomes A n Q° in (z;7) with 

Ql., n (z; 7) = 9°(z; 7)"^ ■ • ■ ^°(*; 7) = Hy)X ■ • ■ < ■ (49) 



To summarize, any complex distribution with support at a single point, eq. (f44|) , admits a 
positive representation, namely, 



p(z) = E -f&lQLJ*) + /^WCdVC*) . (50) 



n=l n[ 



with /3 n given by eq. (f|l]), and 7 given by eq. 

As an illustration we can consider again the distribution of eq. (19). In this case we find 
7 = (4|a| 2 ) -1 and 77 = (5\ = 2\a\ 2 , and thus 





z\ 


2 


4 


\a 


2 



p(z) = \z-2a\ 2 exp(-^). (51) 



As a final application of the results of this section, we can consider periodic probabilities. 
Such probabilities correspond to variables effectively defined in a compact domain and find 
application in the context of compact gauge theories on the lattice. They satisfy, P{x) = 
P(x — na) with (na)i = n^ai where n G Z D is arbitrary and a G is characteristic of P(x). 
Without loss of generality, we may choose a« = 2n. These probabilities do not belong to the 
class previously considered. The normalization as well as the expectation values should be 
taken on a lattice cell {x, < Xj < 2ir, i = 1, . . . , D}. The test functions should be periodic 
and the concept of representation should be modified accordingly: p(z) is periodic on the 
real axis, x is to be integrated on the periodic cell and y on R D . Also instead of equality of 
expectation values of polynomials we demand {exp(irijXj)) p = (exp(injZj)) p for any integers 
rij, j — 1, . . . , D. Assume now that the periodic distribution is a function of the form 

P(x)= ]T P (x-27m), (52) 

where the series is uniformly convergent. Let Pq(z) be a function which is a positive repre- 
sentation of Pq(x) not only on polynomials but also on exponential test functions, and such 
that 
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p( z ) = Po{z - 27m) . (53) 

is uniformly convergent. Then, p(z) is a positive representation of P( readily shown. 

In particular, Pq(x) may be a Gaussian times polynomial and po(z) its positive represen- 
tation found above, since these functions are sufficiently convergent at infinity. Therefore 
the construction given above provides a positive representation for this case too. Another 
example is the periodic version of the one dimensional Gaussian times cosine considered 
above after eq. (Pp: 

P(x) = cos(a;) 2J ex P 
neZ 

p{z) = (6(y - T) + 6(y + T)) £ exp (- (g ~f W H . (54) 

neZ V 21 / 

This example is interesting since it is similar to simplified probabilities considered in the 
literature HT2[ , |15] to model the SU(2) gauge theory in the presence of a Wilson loop, for 
which the complex Langevin algorithm did not work. 



(x — 27rn) s 
2T 



IV. CONCLUDING REMARKS 



We have studied the problem of representation of complex distributions by distributions 
on the analytically extended complex plane. The positive representation problem is of 
immediate interest in some areas of physics: field theory and statistical mechanics. On the 
other hand it also seems a new and interesting field from the mathematical point of view. One 
could consider extending the particular class of complex distributions studied here, namely, 
Fourier transforms of regular distributions analytical at the origin, by allowing as well for 
adding non regular distributions with support outside the origin. Perhaps more interesting, 
and in the opposite direction, one could extend the set of test functions in the definition of 
representation beyond polynomials to insure, for instance, that each probability on C D is 
at most the representation of one probability on K D . From the viewpoint of applications it 
would also be interesting to extend the concept of representations to distributions defined 
on group manifolds since they appear naturally in lattice gauge theories. Our discussion on 
periodic distributions corresponds in fact to the manifold of the direct product of D U(l) 
factors. 
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APPENDIX A: 



In this appendix we will show that (j> n (z) defined in eq. (37) is a null distribution. To 
this end let us introduce the polynomials 
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Q,...,.,j,..jJ^ = 9(2)-%, ■ ■ -d in d* h . ..d* m g(z) . (Al) 
They generalize Qi x ...i n {z) and satisfy the relation 

Qil...inm...jm( Z ) — Q*j 1 ...j m ;i 1 ...i n ( Z ) ■ (A2) 

To prove eq. (|37D, we will use the following Wick theorem: 

Qh...i n (z)Q* h ... jm (z) = £ Q il ... in -,n... jm {z) ■ (A3) 

[ii...i n ;h—jm] 

where the sum is over all possible sets of contractions of the indices i\ . . .i n with the indices 
ji ■ ■ ■ jm- The contraction of two indices i, j gives a factor jSij and removes them from the 
list, e.g., 

Qiii 2 ( z )Q*j( z ) = Qhi 2 ;j( z ) + iKjQiA z ) + iKjQh{ z ) ■ ( A4 ) 

In general there are n\m\/k\(n — k)\(m — k)\ terms with k contractions. Let us apply the 
Wick theorem to Qi 1 ...i n (z)Q* i j n (z)g(z) . Whenever two indices i,j are not contracted we 
will have Qi...;j...{z)g(z) which contains d* and hence is a null distribution. Therefore only 
the terms with all indices contracted contribute and the non null part is 

7" E S iljpl ---5 injpn g(z), (A5) 

pes„ 

where the sum runs over all permutations. After contracting the indices we obtain eq. (p7|). 
K n (D) is the number of ways of choosing n objects out of D allowing repetitions. 
The Wick theorem can be proven by induction. Defining the operator 

Vi =g-\z)d ig {z) = di + Ui, (A6) 

(g(z) is a multiplicative operator here) we have 

Qh...i n (z) =V il ---V ln Q {z), 
Qn...i n ;n... jm (z) = V n ■ ■ ■ V in V* h ■ ■■V) n Q^z) , (A7) 

where Qo(z) = 1. Trivially, [di,V*] = —jSij, thus 

m 

^,..i m W = -E7^g;,.,;.., ro (.), (A8) 

k=l 

where the hat means that the index has been removed from the list. On the other hand 
T>(AB) = (T>A)B — AdB. The Wick theorem holds for n = m = 0. Assuming it has been 
proven up to some (n, m), 

Qh...i n+1 ( z )Q* jl ... jm ( z ) = (Pi n+1 Qh...i n ( z ))Q* jl ... jrn (z) 

= PwiW^W^W) - Qn---in( z )d in+1 Q* jl ... j Jz) . (A9) 
Using that the theorem holds for (n, m) and eq. (|A8|) . 



12 



Qii...i„+i( z )Qj 1 ...j m ( z ) — E Qh-in+V,jl-jm( Z ) 

[ii—i n ;ji—jm] 
m 

+ E E i^ n+ ii k Q ll .., nm .j k .., m ( z ) • (aio) 

[ii...j„;ji...j fc -jm] 

The first term contains all the contractions not involving the index i n+ i, and the second one 
all the contractions involving the index i n+ i, hence the theorem is proven for (n + l,m). It 
is worth noticing that the reverse expansion also holds, i.e., 

Qn..,,,j,..jM = E Qn...i n ( z )Q* jl ... j J z )- (AH) 
[ii—i n ;ji—jm] 

where the contraction of ij now is —7^. 
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